Compression and stretching of a self-avoiding chain in cyhndrical nanopores 
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Force-induced deformations of a self-avoiding chain confined inside a cylindrical cavity, with di- 
ameter D, are probed using molecular dynamics simulations, scaling analysis, and analytical calcu- 
lations. We obtain and confirm a simple scaling relation — / • D ~ 7?"^^* in the strong-compression 
regime, while for weak deformations we find f ■ D = —A{R/Ro) + B{R/ Ro)~^ , where A and B are 
constants, / the external force, and R the chain extension (with Rq its unperturbed value). For a 
strong stretch, we present a universal, analytical force-extension relation. Our results can be used 
to analyze the behavior of biomolecules in confinement. 

PACS numbers: 



Technological advances have made it possible to visu- 
alize and manipulate individual polymer molecules that 
are trapped in nanopores [1, 2, 3, 4]. This enables one to 
test theoretical predictions on confined polymers, shed- 
ding new insights into their static and dynamic prop- 
erties. Besides their technological importance, confined 
polymers are relevant in a number of biological processes. 
In the context of proteins, a newly synthesized protein 
exits the ribosome through a narrow cylindrical pore [5] 
assisted, perhaps, by a tensile force [6]. Both chaperone- 
assisted protein folding and protein degradation by pro- 
teosome involve encapsulation of proteins in cylindrical 
pores, in which they experience compression and stretch- 
ing forces [7, 8]. Furthermore, many bacterial species 
are rod-shaped and even filamentous, and their chromo- 
somes are highly compressed inside the cell and yet well- 
organized [9, 10]. Thus, the study of cylindrically con- 
fined polymers is not only of practical importance but 
also is a first step towards understanding several funda- 
mental biological processes. 

Surprisingly, compared to the progress on the exper- 
imental side, and despite its fundamental and practical 
importance, theoretical understanding of force-induced 
chain deformations in a pore remains behind. In the lit- 
erature, analysis of force-compresion-extension measure- 
ments heavily relies on and is often limited to simple 
scaling arguments [1, 2, 3]. 

The main purpose of this Letter is to present a quanti- 
tative picture of the interplay between confinement and 
force-induced chain deformations (from strong compres- 
sion to almost full stretching). To this end, we combine a 
scaling analysis, molecular dynamics (MD) simulations, 
and a systematic theoretical approach, and obtain forcc- 
compression-extension relations in almost the full range 
of chain deformations. Our results thus provide a quanti- 
tative basis for experiments involving biomolecules (e.g.. 
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FIG. 1: Self-avoiding polymer in a nanopore (a). Simulation 
setting for compression (b) and stretching (c); the two chain 
ends are confined to, but can freely move within, the confining 
"piston" walls. By measuring the pressure exerted on the 
piston walls, we obtain force-compression-stretching relations, 
for a given chain size TV and a pore width D. 



dsDNA and chromosomes) in strong confinement as in 
nano- and micro-channels. 

Consider a self-avoiding polymer with N monomers 
in a cylindrical pore with diameter D, as illustrated in 
Fig. 1, with R being the longitudinal end-to-end dis- 
tance the chain. Owing to the simultaneous presence 
of self-avoidence and confinement, the chain resembles a 
linear string of "compression blobs" [11, 12]; inside each 
blob, the effect of confinement is not significant. Because 
of the clear separation between two length scales (> D 
and < D), the conventional Flory approach may fail to 
capture the correct D dependence of relevant quantities 
(see below). The essential features of the blob-scaling 
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approaches [11, 12] can be reproduced by the following 
"renormalized" free energy [13], 



i?2 



{Nlg)D^ 



B 



D{N/gY 
R 



(1) 



where g is the number of monomers inside a compression 
blob of diameter D, i.e., g ~ [D/a)^/^, A and B are 
constants, and a is the monomer size. (Throughout this 
paper, /3 = l/ksT with ksT being the thermal energy; 
unless otherwise indicated, a is the unit of length.) The 
first term in Eq. 1 describes connectivity of blobs, while 
the second term, which can be rewritten as ~ (N/g) x 
blob density, represents the mutual repulsion between the 
blobs, thus ensuring linear ordering. The second term is 
equivalent to correctly assigning ksT to a binary collision 
of two blobs, independent of their size [11, 14]. 

The free energy in Eq. 1 produces not only the ex- 
pected equilibrium chain size Ro - ND-^/^ [15], when 
J- is minimized, but also the correct confinement free en- 
ergy /3.Fconf = PH^ ^Ro,D)^ N/g ~ ND-^/^, which 
is linear in N [11, 12]. The free energy cost for linear 
ordering is correctly counted as ~ ksT per blob (D- in- 
dependent). Also, the resulting effective Hookean spring 
constant of the chain, /ceff = {d^:F /dR"^) - N'^D^^/^ 
is consistent with the previous scaling result [1, 12]. Note 
that, without the renormalization trick in Eq. 1, Flory 
theory (as employed in [3]; see also [13]) will lead to 
a D-independent (thus qualitatively inaccurate) spring 
constant k^g. 

An important consequence of Eq. 1 is that it pre- 
dicts a universal scaling relation for (external) force- 
compression-extension (i?) relation. 
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To test the validity of Eq. 2 (thus Eq. 1 as well), 
we have carried out molecular dynamics simulations 
using ESPResSo [16] for a wide range of parame- 
ters: D = 4, 5,..., 15 and N = 128,256,512 (see 
Fig. 1). In ESPResSo, monomers (or beads) are con- 
nected by finite extensible nonlinear elastic (FENE) 
(spring) bonds with a Weeks-Chandler- Andersen (WCA) 
potential for excluded-volume interactions (monomer- 
monomer as well as monomer- wall) . We confined the two 
end beads of the chain by "piston" walls at the cylinder 
ends (Fig. l(b)&(c)), while allowing the ends to move 
freely in the plane of the piston wall. By varying the 
wall-to-wall distance (thus R), we obtained the longitu- 
dinal force f{R) to keep R fixed. 

Figure 2(a) summarizes our simulation results, where 
all 24 curves have been rescaled according to Eq. 2. 
The data perfectly collapse onto a single master curve 
(the solid line) for R < Rq. Moreover, in the range of 
0.5 < R/Rq < 1.2, Eq. 2 fits the rescaled data very well, 
thus confirming the scaling prediction of Eq. 1 [see the 
lower-left inset in Fig. 2(a)]. However, as the compres- 
sion becomes strong {R/Rq <^ 1), the blobs break into 
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FIG. 2; [Color online; (a) and (b) have the same color scheme] 
Simulation (dotted lines) and scaling (solid lines) of force- 
compression-stretching of a chain in a cylinder, (a) R/Ro vs. 
/ • D/ksT. All 24 simulation curves collapse onto a single 
universal curve for R < 1.2Ro [N = 128 {D = 4, 5, 6, 7), iV = 
256 {D = 4,5, ...,ll),iV = 512 {D = 4, 5, . . . , 15)]. The 
solid line is the scaling result using Eq. 2 in the range of 
0.5 < R/Ro < 1.2; by a global fit to the simulation, we ob- 
tained 2 A = 2.84 ± 0.03 and B = 2.81 ± 0.01 (i.e., 2Ax B). 
The lower-left inset shows the fitted region. The upper-right 
inset is the log- log plot for R < Rq. (b) R/N vs. |/| for force- 
stretching {R > Rq). The inset shows the linear Hookean 
regime for weak force (see text). The sudden departure from 
the linear relation as R/N 1 reflects the non-universal as- 
pect of the bond potentials in simulations. 



smaller ones [see Fig. l(b)J, and the correct form of free 
energy is /3J^ - (i?3/T/)i/(3'^-i), where i? - A^'^ is the 
radius of gyration of a corresponding unconfined chain 
(with v « 3/5 the Flory exponent) and V ~ D^R the 
confining volume [14]. From this, we obtain 
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[see the upper-right inset in Fig. 2(a)]. 

Importantly, Fig. 2 shows the limitations of Eq. 2; 
when the chain is stretched appreciably from its equi- 
librium length (i? > 1.2i?o), the force-extension curves 
(FECs) do not collapse onto a single universal plot. 
This is not unexpected since an additional length scale. 
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namely the tensile blob size ^ (< D) [11], is relevant at 
high /, as indicated in Fig. 1(c). This also signals the 
onset of a D-independent regime at large /, as the inter- 
actions between the chain and the wall are negligible [see 
Fig. 2(b) and Fig. 3]. 

In Fig. 2(b), we show FECs for R > Rq by rescaling 
R by N. Since R/N denotes the extension per chain 
segment, for the same pore diameter D, the force ex- 
tension curves tend to collapse onto one another, as 
expected. In a recent simulation study [17], the effec- 
tive spring constant k^ff was shown to vary from k^g ~ 
l/ND^-'^ to \/ND^-^ for the intermediate chain length 
(up to iV = 2000). This differs from the scaling result 
/ccff ~ l/ND^/"^ [12], which could be reached only when 
N ~ 10^ [17]. Our results are also in good agreement 
with Arnold et al. and we find the linear Hookean col- 
lapse for fceff - l/7Vi:)"-9±0-2 [/ - k^s{R~ Rq)] inset 
of Fig. 2(b)]. The data cannot be fit using the large N 
result, /ceff l/ND^'^ [18]. 

To further elucidate the subtleties of the FECs in the 



presence of confinement, we have also performed ana- 
lytical, self-consistent calculations following the theoret- 
ical methods described in [19, 20]. Our motivation is 
to derive the FEC relation in the large- limit, which 
can be reached experimentally (if not computationally; 
see above). In the continuum limit, the distribution of 
monomers at r(s) can be represented by the following 
Hamiltonian 

pn ^ ^ [ P{s)ds+'^ I I dsds'5{v{s')-v{s)) 
2a Jo 2 7o Jo 
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where s is the contour length (0 < s < L = No), 
[...] = d[...]/ds, and v the excluded volume parameter. 
Confinement effects will be taken into account through 
the boundary condition that the probability of finding 
any monomer on the wall is zero. 



We replace Ti. by the simpler reference Tii given by PHi = ^ r^{s)ds — Pf z{s)ds and choose an optimal 
value of oi - instead of minimizing the free energy as in the Flory approach - so that R obtained using Tii, i.e., 
R = {zl — zq) = (3Laif/3, coincides with R calculated using H to first order in v [19, 20]. Additionally, we use 
the "ground-state dominance" approximation, which is valid in the large A^ limit [11]. The optimal ai satisfies the 
self-consistent equation (SCE): 
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where ao ~ 2.40 is the smallest (first) zero of the zeroth- 
order Bessel function of the first kind Jq{x), and erf (a;) is 
the error function. When combined with R = (3Laif /3, 
the solution of the SCE determines the equilibrium chain 
length along the pore axis. 

For f — 0, the SCE leads to the equilibrium chain size 
Rq = {{zl — ^o)^)y^o ^ij^y^'j which is consistent 
with the previous result [12, 20]. For / > 0, the SCE 
can be solved numerically. A few comments are in order: 
First, the excluded volume parameter v in the SCE can 
be considered as a fitting parameter (see Ref. [19]). The 
value of V can be adjusted to ensure the best global fit 
to simulations or experimental data. Second, our SCE 
does not reflect the kinetic aspects of force-extensions; 
chain reorientation, allowed in our SCE, is kinetically 
suppressed for strong confinement (It is strictly forbid- 
den in our simulations.) Finally, because of the ground- 
state-dominance approximation, the SCE becomes more 
accurate for large L. In the limit L — > oo, chain reori- 
entation becomes irrelevant in that chain alignment in 
the force-direction is ensured for / > ksT/Ro (— > as 
L co). In this case, v dependence enters the FEC only 



via the combination D'^v and, thus, can be absorbed into 
an effective D (see below). Our focus below is thus the 
large L behavior of our SCE. 

In the large L limit, the SCE in Eq. 5 reduces to 
'^1 (a ^ ot) ~ 0.44-^ f^i^jy . This equation implies that 
ai is a function of /D^/s/^i/a. ^ ai{fD^I^ /v^/^). As 
a result, {AR/N) D^/^v^^^ = (/3ai//3) D^/s/^i/a jg 
function of /D^/ay^i/s^ ^^^^.^ AR ^ R - Rq. Thus, 
the FECs collapse onto each other in a JD"^!^ jv^^^- 
(AR/N) D2/3/uV3 plot. The analytic results justify the 
scaling plot in Fig. 2, where the renormalized Dj^ ^ 
In Fig. 3, we plot ARD^/^/v^^^ as a function 
of jD'^l'^ jv^l^ in the limit A^ ^ oo. All the curves for 
different values of D collapse onto each other (as found 
in the simulations) for the reason described above. (Also 
see the inset in Fig. 3 obtained without D rescaling. This 
suggests that a more strongly confined chain behaves as 
a stronger spring.) 

Three different force-extension regimes are identified 
in Fig. 3: linear (green dotted), power-law (red) and lin- 
ear (cyan). The good linear fit to the very narrow lin- 
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FIG. 3: (Color online) Rescaled force-extension relation in 
the limit iV — » oo. Three distinct regimes are identified and 
presented by the dotted lines (green, red, and cyan), which 
are linear, power-law, and linear fits to the force-extension 
relation in their respective regimes. The power-law analysis 
leads to k^a ~ l/D^^^N. 



rameter range can the linearity be recovered. We note, 
however, that this power-law regime represents a much 
wider parameter space than the initial linear regime. The 
large-/ linear regime is where the force-extension be- 
comes D-independent is due to almost full stretching of 
chain by large force. 

Our results suggest that relaxation dynamics of a 
chain confined in a narrow channel will show several 
distinct relaxation rates (oc fcos) with different D depen- 
dence. This explains the limitations of existing scaling 
arguments based on a single-time scale [1, 2, 3]. The 
predictions based on a scaling approach are reached only 
when is very large. Thus simulations or experiments 
on short chains need to be interpreted with caution. 
Our approach can be extended to incorporate other 
non-trivial but important effects such as chain topology 
(e.g., branched polymers), electrostatics, molecular 
crowding, and chain stiffness [22]. Most importantly, 
the force-compression-extension relations presented 
above can be used to understand the elastic response of 
chromosomes confined in microfluidic devices, as well as 
the fate of proteins in nanopores. 



ear regime indicates a weak D-dependence of the force- 
extension relation [21]. For the power-law fit in the sec- 
ond regime, we have chosen RD^^^ ~ const. -I- {^ fD^/^^ 
with a « 1.5. From this, fceff has been estimated to vary 
as « D~^/^. The seemingly-perfect fit to the 

force-extension in this regime justifies the scaling result 
fceff However, the force-extension in this 

regime is nonlinear, in contrast to what one may expect 
from the scaling approach [12]. Only over a narrow pa- 
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